수치 해석: 스칼라에서 벡터로: 비선형 시스템의 도전

단일 방정식 $f(x)=0$에서 다변수 시스템으로 전환하는 것은 궤도 역학부터 토양 구조 분석까지 복잡한 공학 문제를 해결하는 열쇠입니다. 이제 우리는 직선 위의 단순한 영점이 아니라, $n$차원 공간 내에서 $n$개의 초표면이 동시에 교차하는 지점을 찾는 것입니다.

1. 수학적 구조

비선형 시스템은 각 구성 함수가 미지 벡터 $\mathbf{x} = (x_1, x_2, \dots, x_n)^t$에 의존하는 방정식 집합으로 표현됩니다:

$$f_1(x_1, x_2, \dots, x_n) = 0,$$ $$f_2(x_1, x_2, \dots, x_n) = 0,$$ $$\vdots$$ $$f_n(x_1, x_2, \dots, x_n) = 0,$$

이를 벡터 형태인 핵심 공식으로 요약합니다:

$$\mathbf{F}(\mathbf{x}) = \mathbf{0}$$

여기서 $\mathbf{F} = (f_1, f_2, \dots, f_n)^t$입니다. 개별 함수 $f_i$는 $\mathbf{F}$의 좌표 함수 함수입니다.

2. 해석적 기초 및 연속성

이러한 시스템을 수치적으로 해결하기 위해서는 사상이 잘 정의되어야 합니다. 정의 10.1~10.3은 $\mathbb{R}^n$에서 극한과 연속성이 각 성분별로 결정됨임을 확인합니다.

정의 10.3

함수 $\mathbf{F}$가 $D \subset \mathbb{R}^n$에서 $\mathbb{R}^n$으로 정의된다고 하겠습니다. $\lim_{\mathbf{x} \to \mathbf{x}_0} \mathbf{F}(\mathbf{x}) = \mathbf{L} = (L_1, L_2, \dots, L_n)^t$이 되는 것은 다음 조건이 만족될 때만 가능합니다:

$$\lim_{\mathbf{x} \to \mathbf{x}_0} f_i(\mathbf{x}) = L_i$$ 각 $i=1, \dots, n$에 대해 성립합니다.

$\epsilon-\delta$ 정의를 사용하면, 모든 $\epsilon > 0$에 대해 $\|\mathbf{F}(\mathbf{x}) - \mathbf{L}\| < \epsilon$가 되도록 하는 $\delta > 0$가 존재하며, 이는 $0 < \|\mathbf{x} - \mathbf{x}_0\| < \delta$일 때 성립합니다.

잘못하기 쉬운 점: 노름 독립성

중요한 세부사항: 다양한 노름($\ell_1, \ell_2, \ell_\infty$)을 사용할 수 있지만, 연속성은 특정 선택에 관계없이 동일하다극한의 존재는 $\mathbb{R}^n$ 내 임의의 벡터 노름에 대해 불변입니다.

3. 이론적 복습

정리 1.6: $\mathbb{R}$에서 $\mathbb{R}$로의 함수에 대해서는 연속성을 미분 가능성으로 보일 수 있습니다. 다변수 경우에서는 좌표 함수의 편미분이 존재하고 유계일 경우 연속성이 보장되며, 이는 반복 해법의 전제 조건입니다.

대표 예제: 예제 1

토양 위의 원형 판 문제를 고려해 봅시다. $3 \times 3$ 비선형 시스템을 표준형 $\mathbf{F}(\mathbf{x}) = \mathbf{0}$으로 표현하세요:

$3x_1 - \cos(x_2 x_3) - \frac{1}{2} = 0$
$x_1^2 - 81(x_2 + 0.1)^2 + \sin x_3 + 1.06 = 0$
$e^{-x_1 x_2} + 20x_3 + \frac{10\pi - 3}{3} = 0$

여기서 $\mathbf{x} = (x_1, x_2, x_3)^t$이고 $\mathbf{F}(\mathbf{x}) = (f_1(\mathbf{x}), f_2(\mathbf{x}), f_3(\mathbf{x}))^t$입니다.

질문 1

벡터 사상 $\mathbf{F}(x_1, x_2, x_3) = (x_1 + 2x_3, x_1 \cos x_2, x_2^2 + x_3)^t$에 대해, 왜 $\mathbb{R}^3$의 모든 점에서 $\mathbf{F}$는 연속적인가요?

모든 선형 사상이 연속적이기 때문에 선형 사상이기 때문입니다.

각 좌표 함수 $f_i$가 기본 연속 함수(다항식과 삼각함수)의 합 또는 곱이기 때문입니다.

자코비안 행렬의 행렬식이 일정하기 때문입니다.

모든 $\mathbb{R}^3$를 닫힌 유계 부분집합으로 매핑하기 때문입니다.

질문 2

다음 함수 중 $\mathbb{R}^2 \to \mathbb{R}^2$에서 $(1, 0)$를 제외하고 어디서나 연속적인 함수는 무엇입니까?

$\mathbf{F}(x_1, x_2) = (x_1 - 1, x_2)^2$

$\mathbf{F}(x_1, x_2) = (\sin(x_1 - 1), \cos x_2)$

$\mathbf{F}(x_1, x_2) = (1 / (x_1 - 1), x_2 / (x_1^2 + x_2^2))$

$\mathbf{F}(x_1, x_2) = (e^{x^1}, e^{x^2})$

질문 3

비선형 시스템의 연속성 평가 시 벡터 노름을 유클리드 노름(ℓ₂)에서 최대 노름(ℓ∞)으로 변경하면, 벡터 수열의 수렴 성질에 어떤 영향을 미치나요?

최대 노름이 덜 민감하기 때문에 수열이 발산할 수 있습니다.

유한 차원 공간 내 모든 벡터 노름이 동등하기 때문에 수렴성이 유지됩니다.

극한 값 $L$은 노름 스케일링에 따라 달라집니다.

수열은 ℓ₁ 노름에서 ℓ₂ 노름보다 더 빨리 수렴합니다.

질문 4

뉴턴 방법은 일반적으로 비선형 시스템에 사용됩니다. 그러나 비특이행렬 $A$를 가진 선형 시스템 $Ax = b$에 적용하면 결과는 무엇입니까?

방법이 실패하는 이유는 선형 시스템에 자코비안이 정의되지 않기 때문입니다.

방법은 가우스 소거법으로 축소되며, $n$번의 반복이 필요합니다.

어떤 초기 추정값에서든 정확한 해를 정확히 한 번의 반복만에 도달합니다.

방법은 $x_0$가 $b/n$이 아니면 진동하며 결코 수렴하지 않습니다.

질문 5

함수 $\mathbf{F}(x)$의 극한을 $x$가 $x_0$에 접근할 때 정의할 때, 조건 '0 < ‖x - x₀‖ < δ'는 무엇을 의미합니까?

함수를 정확히 $x_0$에서 평가하고 있습니다.

구멍 난 근처 환경을 고려하고 있으며, 이는 $x_0$ 자체가 극한 계산에서 제외됨을 의미합니다.

원점에서의 거리는 $δ$보다 작아야 합니다.

극한이 존재하려면 시스템이 선형이어야 합니다.